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ABSTRACT 

The effect of primordial magnetic fields on X-ray or S-Z galaxy cluster survey is 
investigated. After recombination, the primordial magnetic fields generate additional 
density fluctuations. Such density fluctuations enhance the number of galaxy clusters. 
Taking into account the density fluctuations generated by primordial magnetic fields, 
we calculate the number of galaxy clusters based on the Press-Schechter formalism. 
Comparing with the results of Chandra X-ray galaxy cluster survey, we found that the 
existence of primordial magnetic fields with amplitude larger than f nGuass would 
be inconsistent. Moreover, we show that S-Z cluster surveys also have a sensitivity to 
constrain primordial magnetic fields. Especially SPT S-Z cluster survey has a potential 
to constrain the primordial magnetic fields with several nano Gauss. 



1 INTRODUCTION 



The origin of large-scale magnetic fields observed in galaxies and galaxy clusters still remains unclear. The most widely accepted 
theory is the astrophysical scenario in which seed magnetic fields are generated by the batte ry mechanism in astrophysical 
phen omena and amplified by the dynamo mechanism in the interstellar or intergalactic medium ferandenburg fc Subramanian 



2005). However, there is uncertainty about the efficiency of the dynamo mechanism IjKulsrud et al.l 



1997 



Giovannini 



2004) 



The recent studies on Faraday r otation measurements of high redsh i ft qu asars suggest the existence of the /^Gauss magnetic 

The existence of such magnetic fields may also 



2008; Bcrnct et al 



2008). 



fields in high-redshift galaxies (jKronberg et al 
challenge the dynamo scenario. 

Anot her candidate for the origin of the galactic magnetic fields is primordial m agnetic fields which is generated in the early 
universe (IWidrowlbool ; iTakahashi et al.ll2005l : llchiki et alJbood : lGiovanninill2008l ) . If primordial magnetic fields existed, Big 
Bang Nucleosynthesis (BBN) and Cosmic Microwave Background (CMB) anisotropy suffer the effect of primordial magnetic 
fields. Therefore, the constraint on the primordial magnetic fields through BBN and CMB anisotropy has been studied by 
many authors. 

After the recombination epoch , primordial magnetic fields generate additional density fluctuations by the Lorentz force 
I Wasserman 19781 ; Kim et al. 19961 ) and many authors have studied their effects on the evolution of the large scale struc 



tures; the redshift-space matter power spectrum, the e p och of reionization, and 21 cm fluctuations IIGopal fc Sethi 



2003 



Sethi fc Subramanian 2005 ; Tashiro fc Sueivama 2006a b ; Schleicher et al. 20091 ; Sethi fc Subramanian! 20091 : Shaw fc Lewis 



2010h . It was found that magnetic fields as small as a few nano Gauss can give strong cosmological impacts. Therefore detailed 



observations planned in near future have the potential to set further constraints on primordial magnetic fields. 

In this paper, we study the effect of primordial magnetic fields on the mass function of galaxy clusters. The additional 
density fluctuations generated by the primordial magnetic fields enhance the formation of galaxy clusters. Thus, number count 
of clusters can constrain their amplitude as well as the standard cosmological parameters such as the amplitude of density 
fluctuations (as) and energy density of matter (fi m ). Specifically we investigate the potential of the cluster number count by 
X-ray observations and Sunyaev-Zel'divich (S-Z) survey. 

We can find massive clusters by the observation of the X-ray emitted from the hot intracluster gas. X-ray-flux-selected 
cluster samples with the calibration between the X-ray temperature and the cluster mass give the mass function of galaxy 
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clusters. Vikhlinin ct al. and iMantz et aJj l |20ich applied it to give a constraint on the cosmological parameters such 
as the equation of state of dark energy. 

The S-Z effect is the scattering of CMB photons by the hot intracluster electron gas (jSunvaev fc Zeldovichlll972l ) and is 
also a powerful tool for detecting galaxy clusters at high redshifts. Combining the S-Z galaxy survey with X-ray or optical 
observations, we can obtain the mass- or redshift-abundance of the number of galaxy clusters. There are many observation 
projects carried out or planned. In particular, Planck and the South Pole Telescope (SPT) is expected to give large catalogs 
of S-Z galaxy clusters and the significant constraints on cosmological parameters. 

The rest of the paper is organized as follows. In Sec. II, we give a description of the density fluctuation generation by 
primordial magnetic fields after the recombination epoch. In Sec. Ill, we study the effect of primordial magnetic fields on the 
mass function by using the Press-Schechter formalism calibrated to numerical simulations. We also compare our result with 
the mass function derived from Chandra observations and obtain a constraint on the strength of primordial magnetic fields. 
In Sec. IV, we calculate the number count for S-Z galaxy cluster surveys and discuss the potential of Planck and SPT to give 
constraints on the primordial magnetic fields. We conclude in Sec. V. Through this paper, we assume a ACDM cosmological 
model with h = 0.72, ft m = 0.3, fi b = 0.05 and a$ = 0.75. 



2 DENSITY FLUCTUATIONS DUE TO PRIMORDIAL MAGNETIC FIELDS 



In this section, we calculate the density fluctuations produced by primordial magnetic fields. First, we make some assumptions 
about primordial magnetic fields. Because our interest is in relatively large length scales, we can assume that the back-reaction 
of the fluid velocity to magnetic fields is small. Therefore, we consider the case where primordial magnetic fields are frozen in 
cosmic baryon fluids, 

Bo (as) 



B(t,x) 



a»(t) 



(1) 



Here Bq(x) is the comoving strength of magnetic fields and a(t) is the scale factor which is normalized as a(to) = 1 at the 
present time, to. For simplicity, we assume that primordial magnetic fields are statistically homogeneous and isotropic and 
have the power-law spectrum with the power-law index n, 

kuk 2 j ' 



<B 0! (fei)B * 3 (fe2)) = ^-5(ki - fa) 



hi 



Ps(k), 



P B (k) 



. /i ll 



(2) 



where ( ) denotes the ensemble average, Bcn(k) are Fourier components of Boi(x), k n is the wave number of an arbitrary 
normalized scale and B n is the magnetic field strength at k n . 

Our interest is to constrain the magnetic field strength on a certain scale in the real space. Therefore, we have to convolve 
the power spectr um with a Gaussia n filter transformation of a comoving radius A, in order to get the magnetic field strength 
in the real space ( Mack et al. 2002T ) . 



Bi 



{B Qi (x)B 0i (x))\ x = 



(2tt) 3 

Substituting Eq. ([2} to Eq. (|3]), we can associate B\ with B 

B 2 



<fkBl [ A 



exp 



A 2 fc< 



(fe n A)-"r((n + 3)/2). 



(3) 



(4) 



(2tt) 2 A 3 

We take h~ Mpc as A throughout our paper. 

Before the recombination epoch, primordial magne tic fields have a damping scale due to the dissipation of the magnetic 



Jjerore tne recombination epocn, primordial magnetic nclds nave a damping scale due to tnc dissipation 01 tne magnetic 
fields by the interaction between electrons and photons ( Jedamzik et al. 19981 ; Subramanian fc Barrow 1 19981 ). As a result, the 



magnetic field power spectrum has a sharp cutoff around the damping scale. The damping scale l/k c after the recombination 
epoch can be related to the magnetic field strength B\ in the power-law magnetic field case as 



fcc = 



143 



VlnGy \0.7 



1/2 



h 2 Q, h 
0.021 



1/2' 



2/n + 5 



Mpc" 



(5) 



in the matter dominated epoch. 

Primordial magnetic fields affe ct the motion of ionized baryon by the Lorentz force even after the recombination epoch 
( Wassermanlll97a : iKim et al.lll996T ). Although the residual ionized baryon rate to total baryon is small after recombination, 
the interaction between ionized and neutral baryon is strong in redshifts considered here. Using the MHD approximation to 
the baryon fluid, we can write the evolution equations of density fluctuations with primordial magnetic fields as, 



(6) 
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S(t,x) = 



V- ((V x B (x)) x B (x)) 
4np b0 a 3 (t) 



d 5 d ad8 d 

-t^tt = -2 — 5- + 4-KG(p h 5 h + p d 8 d ) 
at- 1 a at 



(7) 



(8) 



where pb and pd are the baryon density and the dark matter density, and <5b and 5dm are the density contrast of baryon and 
dark matter, respectively. Solving Eqs. (0 and ((H), we can obtain the power spectrum of the density fluctuations. With the 
assumption that there is no correlation between primordial magnetic fields and primordial density fluctuations, the density 
matter power spectrum can be separated into two parts as 

P(k 1 t) = Pp(k 1 t) + P M (k,t), (9) 

where the first term Pp(k) is originated from the primordial density fluctuations, whose growth rate is proportional to t 2 ^ 3 
in the matter dominated epoch. The second term represents the power spectrum of the density fluctuations produced by the 
primordial magnetic fields. The power spectrum Pyi(k) is written as 

Pu(k)=(^-X( - - tf — ) 2 D M (tfl 2 (k), (fO) 



fi m / \47rft b pcoa 3 (ii) 

where p c o is the critical density at the present epoch, ti is the initial time which is the recombinati on epoch in our calculation , 
DM(t) is the growth rate and asymptotically proportional to t 2 ^ 3 in the matter dominated epoch ( Tashiro fc Sueivamal 2011 ). 
and 



I 2 (k) = (iV ■ (V x Bo(x)) x B (a;)| 2 }. 



(11) 



Under the assump tion of the isotrop i c Gaussian stati stics for primordial magnetic fields, we can rewrite the nonlinear convo- 
lution Eq. (PJ as |Wassermanlll978l : lKim et alJll996h 



I 2 (k) 
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(47r) 2 |fc-fci| 2 

where p is p = fc ■ fci/|fc||fei|. Note that the range of integration of ki in Eq. (|f 2[) depends on k because we assume that the 
power spectrum has a sharp cutoff below l/k c so that ki < k c and |fe — Asi| < k c must be satisfied. 



Eq. (|12p can be estimated analytically in the limit of k/k c <C f as I 2 (k) ~ aB 



2n + 10£.2n + 7 



+ PB^k where a and /3 are 



coefficients which depend on n. Here we employ the fact that the damping scale k c is proportional to B^ as is shown in 
Eq. ([5]). The former term dominates if n < —1.5, while the latter dominates for n > —1.5. Accordingly, if magnetic fields 
have a power-law index smaller than —1.5, the power-law index of density fluctuations depends on that of magnetic fields. 
However, if the power-law index of the magnetic fields is larger than —1.5, the power-law index of the density fluctuations is 
about 4 and the amplitude is decided by their damping scale. 

We introduce an important scale for the evolution of density perturbations, i.e., magnetic Jeans length. Below this scale, 
the magnetic pressure gradients, which we do not take into account in Eq. ([SJ, cou nteract the gravit ational force and prevent 



further evolution of density fluctuations. The magnetic Jeans scale is evaluated as (|Kim et al 



19961 ) 



13. 



VlnG J 



h 2 n m 

0.18 



1/2' 



2/r. 



Mpc" 



(13) 



For simplicity, we assume that the density fluctuations do not grow below the scale, although the density fluctuations below 
the scale are, in fact, oscillating like the baryon oscillation. 



3 MASS FUNCTION AND X-RAY OBSERVATION 

The additional density fluctuations produced by primordial magnetic fields enhance the number of dark matter h alos. In order 



to est imate this enhancement, we use the mass function which is calibrated to fit the numerical simulation bv lTinker et al 
( 20081 ). 
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l(z) = do (1 + z)~ 



b(z) = b (l + z)- a , loga(A) = 



0.75 



-V 



(14) 



(15) 



>g(A/75)y 

where a is the smoothed variance of the density fluctuations with a top-hat window function, Aq, ao, bo and c depend on A 
and A is the overdensity contrast within a sphere of radius _Ra which is related to the halo mass Ma by 



M A = n m ^-R 3 APc A. 

For example, in the case of the halo virial mass, A is A = 178 in the matter dominated epoch. 



(16) 
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Figure 1. The integrated number counts of X-ray galaxy clusters. The left panel shows the number count for the low redshift bin 
z = 0.025 — 0.25, while the right panels shows that for the high redshift bins z = 0.35 — 0.90. In both panels, the solid, the dashed and 
the dotted lines repr esent the numbe r coun ts for the case of B\ = 10 nG, Bx = 8 nG and B\ = nG, respectively. We also plot the 
results of Chandra in Vikhlinin ct al. For comparison, we put the number counts for the B\ = nG case with erg = 0.85 as the 

thin solid line. 



Our interest is to examine the potential of the X-ray galaxy cluster observation to give constraints on primordial magnetic 
fields through the halo m ass function. In t he X -ray observation, the halo mass for high A is more robust than for low 
A. Therefore, according to Vikhlinin et al. ( 20091 ') . we set A = 500. The observational luminosity threshold gives the mass 
threshold for observed halos. Therefore, the number count of halos over the luminosity threshold corresponds to the integrated 
mass function, 

N > M * = L™m> (17) 

where Ml is the mass threshold corresponding to the luminosity threshold. 

In Fig. [I] we plot the integrated number count as a function of the mass threshold for each primordial magnetic field 
strength with the magnetic field spectral index n = —2.8. The left panel in Fig. [T] is for the low redshift bin, z — 0.025-0.25, 
while the right panel is for the high redshift bin, z = 0.35-0.9. Although the additional density power spectrum induced 
by primordial magnetic fields dominate the primordial power spectrum on smaller scales than 1 Mpc, the additional power 
spectrum can enhances erg to 0.8 for B\ = 10 nG. For reference, we put the integrated number count for the ACDM model with 
B\ = nG and erg = 0.85. The existence of the primordial magnetic fields lift up the mass function on small scales, because 
the additional density fluctuations produced by primordial magnetic fields have a blue spectrum. As a result, compared with 
the case with B\ = nG and erg = 0.85, while the number count in the case with primordial magnetic fields is small on large 
scal es, it is more enhanced on small scales. 

Vikhlinin et~ai] J2009) obtained galaxy cluster mass functions at two redshift range using Chandra observation data and 



concluded that the obtained mass functions are in good agreement with the cosmological model with erg = 0.75. We compare 
the theoretical mass functions with the data with the error bars due to the Poisson uncertainties in Fig. [T] From this figure, 
we can conclude that Chandra observations rule out the primordial magnetic field strength, Bx > 8 nG at roughly one-sigma. 



4 S-Z NUMBER COUNTS 



The S-Z effect is caused by the scattering of CMB photons with electrons in hot ga s in galaxy clusters. The change of the 
CMB intensity with the frequency v by the S-Z effect is expressed, in the R-J limit, as jSunvaev fc Zeldovich|[l972l ; iBirkinshawl 
19991 ) 



I V (B) = 2v 2 TcMng{x)y{e), (18) 

where Tcmb is the CMB temperature and g(x) is the S-Z effect spectral shape given by g(x) = s 2 e a; (a;/tanh(a;/2)-4)/(e :,; -l) 2 
with x — 2-kv/Tcmb- The Compton y-parameter is given by the integral of the electron gas pressure along the line of sight at 
9 

y{6) = [dl f g ^^n e a T , (19) 
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where T e is the electron gas temperature, n e is the electron density, <tt is the Thomson scattering cross-section, m e is the 
electron mass, /i e is the mean mass per electron ri e = 1.143 and f g is a gas fraction in a galaxy cluster which we set f g = 0.12 
i Mohr et all 1 19991 ). In the S-Z cluster survey, it is assumed that the S-Z cluster is a point like source within the telescope 
beam. Therefore, we consider the total flux density from a cluster at redshift z by integrating the cluster surface, 



2 Y 
2v TcMBg(x) — , 



S v =2v*TcMBg(x)i^, (20) 
where D a is the angular diameter distance to the cluster at z and Y is the integrated {/-parameter over the cluster surface, 



Y 



dtt y(9). 



(21) 



In order to calculate Eq. (|21[) . we need the e lectron density profile in a galax y cluster. We take the assumption that the 
electron density profile is isothermal /3-profile ( Cavaliere fc Fusco-Femianoll 19761 ) . 

M 



n e (r) 



no 



H2 , r < R v 

r > i? v 



(22) 



where R c is the core radius of galax y cluster which is related to the virial radius R v with the parameter s(z) as R c = R v /s(z) 
Following iKomatsu fc Seliakl (|2002r i, we set 

-0.2 



s(z) 



10 

1 + 2 



M 



M*(0) 



(23) 



where M*(0) is a solution to c(M) = <5 C at the redshift z — 0, where <5 C is the critical density contrast for collapsing. 

Taking the assumption that the galaxy clusters are spherical and in the hydrodynamical equilibrium, we can relate the 
virial mass M v with the virial radius, 



M v = —pM(z)A v (z)Rt, 



(24) 

where A v (z) is the overdensity contrast for virialization jNakamura fc Sutol 1997 ). 

We introduce the electron density weighted average temperature (T e ) n , which is defined by (T e ) n = J dln e T e / J dln e . 



Battve fc W cllcr (2003|) has obtained (T e ) n under the isothermal assumption in the ACDM model, 



<Te)„ = T* keV 



Mvir 



-i 2/3 



lO^/i-iM 
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2T 1/3 
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, (i-n m (z)) 

' A c ( 2 ) 



(25) 



where T+ is the tempera ture normalization factor a nd they adopt = 1.6 for agreement with numerical simulation works in 



Bryan fc Normanl || 19981 ) and IPierpaoli et all (|200ll ). 



Using the /3-profile assumption with the electron density weighted average temperature given by in Eq. (|25|l , we can write 
the y-parameter as 



where £(9) is the projected profile of the electron density, 



/Lt e m p 



.1 tan 



-a T ao), 



s 2 -0 2 /0* 

T+WJef 



1/2 



tan 



(26) 



(27) 



In actual observation s, the finite be am size of telescopes causes the beam-smearing effect. This effect can be accounted 
by modifying Eq. (J2TJ) to (|Bartlettll200ol ) 



Y 



dQ y(8)B(8). 



exp[ 



(28) 

7(2a 2 )] with a b — 0fwhm /V8 In 2 



Here we assume that the beam profile is described in a Gaussian form B(9) 
where Sfwhm is the full- width-half-maximum (FWHM). 

The parameter Y depends on the mass and the redshift of galaxy clusters. Therefore, giving the flux limit Sum of the 
observation, we can obtain the limit mass of galaxy clusters Ml at each redshift. We show Ml for Planck and SPT in Fig. [2] 
We present the parameter value for each observation in Table [1] Planck will cover the full sky and the Planck sensitivity is 
14 mjy at 100 GHz. SPT covers Att = 4000 degree square and the SPT sensitivity is 0.8 mjy at 150 GHz. SPT has a better 
sensitivity than Planck. As a result, Ml for SPT is lower than for Planck. 

The combination between the S-Z galaxy cluster survey and the follow-up optical observation enables us to obtain the 
cluster number count for redshift bins. In Fig. [3] we show the number count of galaxy clusters with mass higher than the 
limiting mass shown in Fig. [2] 
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redshift 

Figure 2. The limiting mass for each redshift. The solid line is for Planck and the dashed line is for SPT. The each experimental 
parameters are in Table. [T] 
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4000 



Table 1. The experimental parameters for S-Z surveys. 

In both panels of Fig. O we set Az — 0.1 and the primordial magnetic field spectral index n = —2.8. 

Fig. [3] shows that the predicted number counts of Planck and SPT are almost same in low redshifts. This is because, 
although Planck has less sensitive to small galaxy clusters than SPT, Planck's full sky survey area increases the number of the 
observed galaxy clusters, and vice versa. However, in high redshift, since the number of large-mass clusters rapidly decreases, 
the number count for Planck become much lower than for SPT. 

Primordial magnetic fields generate additional density fluctuations in small scales and bring the early structure formation. 
Therefore, the difference from ACDM cosmology due to the existence of primordial magnetic fields is emphasized in the SPT 
observation, especially in high redshifts. Even the primordial magnetic fields with B\ = 6 nG amplifies the number count in 
high redshifts by 50 % for the SPT sensitivity, while, for Planck, such primordial magnetic fields cannot bring a significant 
amplification. 

For reference, we put the number counts for the B\ = nG case with as — 0.85 as the thin solid line in both panels in 
Fig. [3] Although as in the case of B\ = 10 nG corresponds to 0.8, the power spectrum in the case of B\ — 10 nG has larger 
amplitudes in high fcs than in the case of as = 0.85 without primordial magnetic fields. This results in the fact the number 
of galaxy clusters with small mass is larger in the case of B\ — 10 nG than in the case of as = 0.85. Therefore, the number 
counts for B\ = 10 nG exceed the ones for as = 0.85 in the low redshifts where both Planck and SPT are sensitive to low 
mass clusters as shown Fig. [2] In particular, the number counts of SPT for B\ — lOnG is almost same as for erg = 0.85 even 
in high redshifts, because SPT has small limiting mass. 



5 CONCLUSION 

In this paper, we have studied the effect of primordial magnetic fields on the galaxy survey by X-ray and S-Z observation. The 
primordial magnetic fields generate additional density fluctuations which has a blue power spectrum. Therefore, the number of 
galaxy clusters, especially small ones, is enhanced. X-ray and S-Z survey can directly observe this enhancement. Nano-Gauss 
primordial magnetic fields bring observable enhancement of the number count by the order of factors. 

For X-ray cluster surveys, we have used Chandra's result to put a constraint on the amplitude of primordial magnetic 
fields. We have found that Chandra's result rules out the existence of primordial magnetic fields with B\ > 10 nG at roughly 
one-sigma level. 
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Figure 3. 
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Bx = 6 nG 
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The integrated number counts of S-Z galaxy clusters. The left panel is for Planck, and the right panels is for SPT. In both 
solid, the dashed, the dashed-dotted and the dotted lines represent the number counts for the case of B x = 10 nG, Bx = 8 nG, 
and Bx = nG, respectively. For comparison, we put the number counts for the B x = nG case with as = 0.85 as the thin 



S-Z cluster surveys also have a sensitivity to constrain primordial magnetic fields. Especially the observation like SPT 
which has small limiting mass with 1 arcmin angular resolution is a good probe of primordial magnetic fields. We have found 
that the combination with high redshift optical surveys has the potential to put the constraint on the fields of nano Gauss 
order. 

In this paper, we consider only primordial magnetic fields with n = —2.8. The power spectrum of the density fluctuations 
generated by primordial magnetic fields has a dependence on the spectral index of the primordial magnetic fields. The large 
spectral index induces the large amplification of the density fluctuations on small scales and increases the mass function for 
small-mass clusters. For example, Bx = 4 nG and n = —2.5 amplifies the number count in high redshifts by 50 % for the SPT 
sensitivity, comparing with the number count without primordial magnetic fields. This amplification is same as in the case of 
Bx = 6 nG with n = —2.5. Therefore, SPT has the potential to put the strong constraint on the primordial magnetic fields 
with large n. 

In our calculation, we ignore the effect of primordial magnetic field on the structure of a halo. However, in order to obtain 
a highly accurate constraint on primordial magnetic fields, it is necessary to study the modification on the electron density 
profile and t he re lation between the X-ray temp erature and the cluster mass by primordial magnetic fields. For example, 



ra y terru> 

Zhand (|2004h and lGopal fc Rovchowdhurvl ( 20101 ) pointed out that magnetic fields with several /^Gauss in a halo modify the 



electron density profile and this modification change the S-Z effect signal. The adiabatic contraction in the halo formation 
easily amplifies the order of nano Gauss of primordial magnetic field strength to the order of /i Gauss. Taking into account 
such effects, we will study the constraint on primordial magnetic fields through X-ray and S-Z surveys in the future. 
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